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Abstract. In the rest frame of the nucleus, shadowing is due to hadronic fluctuations of the incoming 
virtual photon, which interact with the nucleons. We expand these fluctuations in a basis of eigenstates of 
the interaction and take only the qq component of the hadronic structure of the photon into account. We 
use a representation in which the qq-pah has a definite transverse size. Starting from the Dirac equation, 
we develop a path integral approach that allows to sum all multiple scattering terms and accounts for 
fluctuations of the transverse size of the pair, as well as for the finite lifetime of the hadronic state. First 
numerical results show that higher order scattering terms have a strong influence on the total cross section 

<Jt ot A - The aim of this paper is to give a detailed derivation of the formula for the total cross section. 

PACS. 11.80.Fv Approximations (eikonal approximation, variational principles, etc.) 11. 80. La Multiple 
scattering 13.60.Hb Total and inclusive cross sections (including deep-inelastic processes) 25.20.Dc Photon 
absorption and scattering 
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1 Introduction 

The experimental observation that the total 7*-nucleus 
cross section at small Bjorkcn-x, Xb, is smaller than A 
times the 7*-nucleon cross section, 



(1) 



is called shadowing. Many theoretical efforts have been 
devoted to understand this phenomenon quantitatively. A 
broad review of the experimental and theoretical situation 
can be found in pj. Depending on the reference frame, 
different physical pictures arise. In the Breit frame, the 
nucleus appears contracted and parton fusion leads to a 
reduction of the parton density at low Bjorken-x 
A very intuitive picture arises in the rest frame of the nu- 
cleus, where shadowing may be understood qualitatively 
in the following way: The virtual photon fluctuates into a 
hadronic state |6|-(ljl that interacts with the nucleus at 
it's surface and the nucleons inside have less chances to 
interact with the photon. Thus, the total cross section is 
smaller than expected naively. 

When we want to describe shadowing in the rest frame 
of the nucleus, we have to choose an appropriate basis 
in which the hadronic fluctuation is expanded. Since the 
hadronic states must have the same quantum numbers as 
the photon, it is reasonable to write the physical photon as 



a superposition of vector mesons. This idea leads to the 
(generalized) vector meson dominance model (G)VMD. 
The fluctuation extends over a distance called coherence 
length 

Lc = Q 2 + m* ' ( 2 ) 

where Mx is the invariant mass of the fluctuation, v is the 
energy of the photon and Q 2 it's virtuality. For small xb, 
this length can become much larger than the nuclear ra- 
dius Ra- In this hadronic basis |lq| , the shadowing correc- 
tion is given by the Karmanov-Kondratyuk-formula [^6| , 
i. e. 



7*A 



A<j. 



4n 

T 1 



(T) 



d 2 b I dM 2 x 



d 2 a(j*N -> XN) 



dM\ dt 



(3) 



Fl(L c ,b), 



t=o 



in the double scattering approximation. For the case, L c ~ 
Ra, the finite size of the nucleus has to be taken into 
account. This is encoded in the nuclear form factor, 



1 



A(T) 



dz n A {b,z)e iz/L ' 
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with 



<T) = - / d 2 bT 2 (b). 



(4) 



(5) 
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Here, T (b) = dz tia (b, z) is the nuclear thickness, 
i. e. the integral of the nuclear density over the direc- 
tion of the incident photon and b is the impact parame- 
ter. However, formula (||) takes only the double scattering 
term into account. When we want to calculate corrections 
from higher order scattering terms, we are faced with the 
problem, that the vector mesons are not eigenstates of 
the interaction and processes like XN — ► X'N, where the 
meson is scattered into another state, are possible. 

This problem may be solved by using the eigenstates 
of the interaction as basis ]r7fl . For very low Bjorken- 
x, xb < 0.001, such an eigenstate is a quark-antiquark 
pair with fixed transverse separation p. The separation 
is frozen during the propagation through the nucleus, be- 
cause of Lorentz time dilatation. At higher values, xb ~ 
0.01 or xb ~ 0.1, the finite size of the nucleus will be im- 
portant and we have to introduce the nuclear formfactor 
Fa (L c ,b). This is a serious problem, because no system- 
atic way is known to implement Fa (L c , b) into higher or- 
der scattering terms than double scattering, cmp. eq. (||). 
Even worse, we do not exactly know what Fa is, because 
the coherence length, eq. (||) , depends on the mass of the 
hadronic state Mx- This quantity is well defined, when 
we use the hadronic basis, but for two quarks with fixed 
transverse separation, no mass is denned. 

A solution to these problems has been proposed in JOl 
and numerical calculations have shown the importance of 
multiple scattering, especially for heavy nuclei. For lead 
it gives a correction about 50%. However, the formula for 
<Jt ot A was given without derivation. Therefore, the aim of 
this paper is to give a complete derivation and to point 
out the approximations. 

Before we start with the derivation, we briefly sketch 
the idea of the approach. When no higher Fock-states are 
taken into account, the total cross section cr] ot A is the cross 
section for production of a gq-pair in the field of the nu- 
cleus. This means, we calculate DIS from the elastic scat- 
tering of the qq component of the virtual photon off the 
nuclear target. The total cross section for pair production 
on a single nucleon may be written in the form || 



Here, Zf is the flavor charge, mt the mass of a quark of 



7*JV 



dX 



<Pp [\$T(ep)\ 2 +\$L {ep)\ 2 )a"{p) 



(6) 

with the total cross section er^ (p) for the scattering of 
the pair off a nucleon. The probability for the virtual pho- 
ton to fluctuate into a qq-pair of transverse separation p 
is described by the transverse and longitudinal light-cone 
wavefunctions, summed over all flavors, colors and spin 
states 



\$t(sp)\ = 

\ 2 



\®l (ep)\ 



^5> 2 {(1-2A(1-A)) £ 2 ^M 
7r| / i 

m\K Q {spf}, 

N f 

} J2zjQ 2 X 2 (l~\) 2 K (e P ) 2 . (8) 
(A) f=\ 



24a em 



(7) 



1/137 and e 2 = A(l- X)Q 2 + m 2 f , A is the 



flavor /, ot e 

light cone momentum fraction carried by the quark. Kq 
and K\ are the MacDonald functions of zeroth and first 
order, respectively. We point out, that the pair is created 
electromagnetically in a color singlett state, but interacts 
with a nucleon via pomeron exchange. 

For a nuclear target, the total cross section may be 
written in eikonal form H 



■y'A 



d\ / d 2 p [\<P T (ep)\ + \$L(ep)\ 



x 2 / d 2 b 



l-expf-^W 



(9) 



if the transverse separation p is frozen, i. e. at very small 
xb- In this approximation, one splits a fast oscillating 
phase factor from the wavefunction ip(r) = e lkz ip(r) and 
obtains for the slowly varying part (p an equation of the 
form 

a ~ 1 ^ (10) 



OZ 



const. 



V 



In the target rest frame, the interaction is given by a color- 
static potential V. However, anticipating that all depen- 
dence of the interaction will be absorbed into the dipole 
cross section cri^(p), we use an abelian potential. Strictly 
speaking, this is justified only in the case of an electron- 
positron pair propagating in a condensed medium, but our 
most important aim is to demonstrate explicitely, how to 
treat fluctuations of the transverse size of the pair, for val- 
ues of Xb , where the transverse size is not yet frozen. Since 
the influence of the potential will be expressed in terms of 
scattering amplitudes and because our final result inter- 
polates between Aa] ot N and eq. (^|), we assume, that our 
results hold also for the case of a nonabelian potential and 
in the presence of inelastic processes. 

The Laplacian acts only on the transverse coordinates 
and is omitted in the eikonal approximation. Then, eq. 
( |To|) is easily integrated. The multiple scattering series is 
summed like in Glauber theory |y|. This is possible, be- 
cause the typical distance between two nucleons inside a 
nucleus is roughly 2 fm, while the gluon correlation length 
is much smaller, presumably ~ 0.3 fm. This approxima- 
tion has been studied for the case of QCD by Mueller [|o| 
with the result that the scatterings off different nucleons 
inside the nucleus are additive. After averaging over the 
medium, one obtains eq. (Q). 

The idea is now to keep this Laplacian, because it 
describes the transverse motion of the particles in the 
pair. Taking into account this motion, we can correctly 
describe the effective mass of the fluctuation in a coor- 
dinate space representation. The phase shift function has 
then to be replaced by the Green function for eq. (^0|). 
We write this Green function as a path integral and aver- 
aging over all scattering centers yields an effective Green 
function W with an absorbtive optical potential V op t — 
—iVqq (p) ha {b, z) /2. All the details will be given in the 
next section. 
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Fig. 1. A cartoon for the shadowing (negative) term in (|ll|). 
The Green function W(p 2 , Z2\p 1 , zi) results from the summa- 
tion over different paths of the qq pair propagation through 
the nucleus. 



To make the influence of multiple scattering as clear as 
possible, we represent the total cross section in the form 



j* A 

n t 



An 1 * 1 * 
Aa tot 



a tot , 



(11) 



similar to eq. (Q), where crj™ t * is an interference term that 
accounts for multiple scattering. This term, illustrated in 
fig. (P, is due to destructive interference between pairs 
created at different coordinates z and contains also the 
contribution of the transverse motion of the pair to the 
coherence length. It is described by the Green function W 
and avoids the appearance of the undefined quantity Mx- 
The Green function also takes the finite size of the nucleus 
into account. 



2 Derivation of the formula 

In order to derive explicit expressions for the terms in eq. 
(|ll|), we start from the general expression for the cross 
section for the production of a gg-pair, 



da^ A = Z]a em \M fl \ 2 S (is - p° q ~ p° q ) 



d 3 p q d 3 p q 
[2n)\2pp Pq > 



Z f a e 



\M fi \ 



d 2 p±,q d 2 p±,qd\ 

(2tt) 4 Ais 2 A (1 - A) 



(12) 
(13) 



where Z 2 is the flavor charge, a em = 1/137 and the matrix 
element is given by 



M fi 



(14) 



We choose the z-axis to lie in direction of propagation of 
the photon. The photon's momentum is denoted by k and 
it's energy by is, p q is the momentum of the quark and p q 
the momentum of the antiquark. In eq. ( |i"3[ ) we have intro- 
duced the energy fraction A = p\jv. In the ultrarelativistic 
case we are considering, pair production takes place pre- 
dominantly in forward direction and therefore we distin- 
guish between the longitudinal direction (z-direction) and 



the transverse directions. The longitudinal momenta are 
large compared to the flavor masses, p z 3> rrif, and to the 
perpendicular momenta, |pj_| 2 ~ m 2 . In our approxima- 
tion, we keep terms of order m,f/p z only in exponentials 
and neglect them otherwise. All higher order terms are 
omitted. 

Note, that the pair is created electromagnetically in a 
color singlet state, therefore we have the factor Zfa em in 
eq. ( |l2| ) , but we describe the with the nucleons an abelian 
potential </> (r) . The particles in the pair move in a poten- 
tial U (r) that is a superposition of the potentials of all 
nucleons, 

A 

U(r) = J2<P(r-r j ). (15) 
i=i 

The vector r,- runs over all positions of the nucleons. 

The wavefunction of the quark fullfills the Dirac equa- 
tions and is an eigenstate with positive energy, while the 
antiquark is represented by an eigenstate with negative 
energy, 

(p° q -U(T)-m f p + icf V)«P,(r) =0, (16) 
{-P°g -U(r)- mfP + la ■ V) 9 9 (r) = 0. (17) 

No interaction between the quark and the antiquark is 
taken into account and therefore, the two equations de- 
couple. The wavefunction of the quark & q (r) contains an 
outgoing plane wave and an outgoing spherical wave in 
it's asymptotic form, while the wavefunction of the anti- 
quark contains an incoming spherical wave and an incom- 
ing plane wave. We transform these equations into sec- 
ond order equations by applying the operator (pp q — U (r) 
+rrif(3 — iot ■ V) on the first and the corresponding oper- 
ator on the second equation, as described in pif . When 
we omit the term quadratic in the potential, we obtain 



A + | Pg | 2 - 2p° q U(r) + ice ■ (VC7(r))J V q (r) = 0, (18) 
A + | Pg -| 2 + 2p° q U(v) + ia ■ (VU(r))) * q (r) = 0. (19) 



The solutions may 
Sommerfeld-Maue [22 

!p,(r) = e tp " r ( 1 



>proximately be written as Furry- 
|23| type wavefunctions, 



V F q (r)u( Pq ,X q ), (20) 



W q (r) = e-^ r (l + H • V) F g (r) v ( Pq , X q ) . (21) 

Here, u (p q , X q ) is the free spinor with positive energy and 
polarization X q . It satisfies {jj> q — m/) u (p q , X q ) = and 
similarly {$qY rn,f)v (p q , X q ) = 0. In Dirac representation 
they read 




(22) 
(23) 



4 



J. Raufeisen et al.: Nuclear Shadowing in DIS at Moderately Small xb 



The three Pauli spin matrices are denoted by <x and the 
Pauli spin state referred to the rest frame of the particle 
is Xqi or Xq respectively. This means explicitely s • er \q = 
XqXq and s-cr \q = ~X q Xq with a spin vector s normalized 
to unity and X q , Xq — ±1. The functions F q and F q have no 
spinor structure any more. They contain all dependence of 
the potential and have to be calculated for a given U (r) . 
They fullfill the equations |l| 

(A + 2ip q ■ V - 2p° q U (r)) F q (r) = 0, (24) 
(A - 2tp q ■ V + 2pp (r)) Fq (r) = 0, (25) 

with boundary conditions F —> 1 for the quark and for 
the antiquark as z — > oo. Note, that it is essential to take 
the correction proportional to a in eq. ( p0| ) and ( pi]) into 
account, although these terms seem to be suppressed by 
a factor of l/p°. It turns out, that when we calculate the 
matrix element of the current operator, (a • e.) u, between 
free spinors, the large part cancels and we are left with a 
contribution of the same order as produced by the correc- 
tion term. Thus, the terms proportional to a may not be 
neglected in the matrix element, although they give small 
corrections to the wave functions. 

In order to remove the dependence on the transverse 
momenta from the phase factors in eq. ( pp| ) and we 
rewrite the solutions in the form 

-r , \ i\n \r ( \ i a ~* a l- 



W q (r) = e 



-iln- Iz / -, i« _ ac± 
- o l|P « |Z ' H • V 



2p? 2p% 



P-L,. 



o viPg I - ) ^ ( r ) u (P?> A ?) . ( 27 ) 



with 



ijj q ( r ) = e* J -«- r - L e- i ^ p « (r) , (28) 

V-g(r) = e - i P- t '*- rx e < fl p *l-"«'*>*F ff (r) , (29) 



and |p | = y ]3q — for the quark and the antiquark 

respectively. In the following, we neglect the terms con- 
taining a z in eq. ( pq ) and (|27j), because they are of order 
O (l/pqPqty. The functions i/) q (r) and ipq( r ) Wln pl &v the 
role of effective wave functions for the quarks. 

The phase factors combine in the matrix element ( |l4|) 
to the minimal longitudinal momentum transfer, 



qr n = k-\ Pq \-\pq 

and we obtain 

M fi = [ d 3 re* q ™" z 



/ 



2v 2p , g 2p ,9 ' 



(30) 



u r (p g , Ag) ( a e + a ■ ^ ^ a ■ e 



iV (r ? ) - p ± g , 
-a • e a ■ — v (p q , X q ) 



2p°q 



x < (r,) ipq (r q ) 



r=r a =r s 



(31) 



Here, the coherence length, L" lax = l/q™ m , comes into 
the game as an oscillating phase factor. However, l™ ax 
does not depend on the transverse momenta. Their in- 
fluence on the cross section is encoded in the rest of the 
wave functions, eq. (|2^) and (27). The operator V (r q ) acts 



only on the variable v q and the operator V (r^) only on 
r q . After the derivatives have been performed, the whole 
integrand has to be evaluated at r — v q = Yg- 

With the representation (|2^ ) , (^3|) we obtain after 
some algebra within the demanded accuracy 



M 



fi 



d 2 r± 



1 



x x\{ m f° ■ £T 



+i (1 - A) a ■ e z e T ■ V (r ±i? ) + iXcr ■ e z e T ■ V (r_i_ jg ) 
+ (1 - A) (e„ x e T ) • V (r x , ? ) — X(e z x e T ) • V (r x , g ) 



+ 2QA(1-A) \xq 



x ip* q (r_L,g, z) ipq (r ±i g, z) 



(32) 



Some details of the calculation can be found in p4| . The 
unit vector in z-direction is denoted by e z . The polarisa- 
tion vector ex corresponds to transverse states of the 7*, 
while the last term in the curly brackets is due to lon- 
gitudinal polarisation. There are two remarkable aspect 
concerning this last equation. First, it does not contain 
any derivative with respect to z any more, because of 
the transverse nature of the polarization vector ex and 
of e z x er- Second, all dependence on the transverse mo- 
menta in the spinor part has cancelled. 

From the eq. (H) and fl25| ) and from the definition 
of tp, eq. ( p8|) and (P9|), we can obtain an equation for 
ip. Assuming that these functions are only slowly varying 
with z, we omit the longitudinal part of the Laplacian and 
obtain 

(A, " 

Wq 



i—^ q (r x , z) =[^ + U (r ± ,z) ) ^q (t ± ,z) . (34) 



Our ansatz yields two dimensional Schrodingcr equations, 
where the z-coordinate plays the role of time and the 
mass is given by the energy. The Laplacian A± acts on 
the transverse coordinates only. The functions i\) q (rj_ , z) 
and %l>q (r j_ , z) become two dimensional plane waves for 
z — > 00, up to a phase factor that cancels in the square 
of the matrix element. It should be mentioned that the 
kinetic energy for the antiquark has a negative sign. This 
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results from the fact that solutions of the Dirac equation 
with negative energy propagate backwards in time. The 
Laplacian in ( |33| ) and (]34|), account for the transverse mo- 
tion of the pair in which we are especially interested. The 
functions t/j q (r±,z) and ipq ( r _i_j-z) in the matrix element 
( |32| ) may now be expressed in terms of the Green functions 
for (133) and (B4J) and its asymptotic behaviour, 



(35) 



ipq (r_L,2, z 2 ) = J d 2 r^s G q (r x , 2 , z 2 | r x ,i, 



ip±, q -r±,i p -i(\p q \-Pz, q )z 



ipq (r± i2 , z 2 ) = J d 2 r ±A Gq (r± t 2,z 2 \ r X) i, Zoc) 

x e -iPX,fPJ.,lg<(|p(f |-P*,5)«aa_ 

Now we use the expression for do 1 A , eq. (|l^), and 
with the matrix element (|3^) and the last two relations, 
we obtain 



da 1 A = Z 2 f a em 



dz 



d Z > e iir n •(*-*') 



efr± d 2 Tq d 2 Tq J d z r' x d 2 r' q d 2 T'- q 
x 0(r± iq ,r± ! q)G*(r± iq ,z\T q ,z 00 ) 

Gq {?±,q, Z | Tq, Zoo) 

r±=r±, q =r±,q 

x 0*{v'^v' ±tq )G q {vl q ,z'\T' q , Zoo ) 

G*q { T> L,qi Z ' I T 'qi Z oo) 



x e 

X 



'{ T q - T i)-P-L,g e i { T 'q- T e)'P-L,5 



d 2 p± iq d 2 p± i qdX 



(2tt) 4 A 2 (l- A)^ 2 
For convenience we have introduced the operator 
0(r ±tq ,r ±tq ) 



= x\{ ™ f cr ■ e T 



+ i (1 - A) a ■ e z e T ■ V (r Xj g) + «Aer • e z e T ■ V (r x>g ) 
+ (1 - A) (e, x e T ) • V (r X)f ) - A (e a x e T ) ■ V (r x , g ) 



(37) 



+ 2QA(1-A) \ X q. 



(38) 



In order to obtain the total cross section, we integrate over 
P±,q an d P±,q- The exponential factors give a <5-function 
that enables us to perform the integrations over all the 
ts. Note, that O does not depend on r. We get 



7* A 



J 7 J 



d\ 



A 2 (1 - A)' 



dz / dz'e^"^ 7 - 7 -') \ d 2 



d 2 r^ 



x 0(r x>s ,r x , ? )C9* (ri )9 ,ri iff ) 
x G g (rj_ !g , z' |r Xjg , z) 



(39) 



Instead of four propagators we are left with only two, be- 
cause we have used the convolution relation 

G(r± t 2,Z2 \ t± i i,zi) 

d 2 r ± G {r ± . 2 ,z 2 | r x , z) G (r x , « | r x> i, z x ) . (40) 



In order to derive an expression that is convenient for 
numerical calculations, we make use of the path-integral 
representation of the propagators in eq. (B9L) . They read 



G(ri,*'|r x ,*) 

=JvTC^hJd^(±^T 2 -U(T,0) } (41) 



T>t exp { i I ±^-t 2 g?£ 



3=1 



(42) 



where the upper sign corresponds to the quark and the 
lower to the antiquark. In this expression, r is a function 
of £. The derivative with respect to £ is denoted by t. 
Obviously, the condition 



G(r' ± 



z \r±,z) 



S (2) (ri - r x ) 



(43) 



has to be fullfilled and further it must obtain r (z) = 
and r (z 1 ) = r' ± . In eq. (|42j), we have introduced the phase 
shift function X(r — r_^jj. The step function (z) is 
for z < and 1 for z > 0. As mentioned before, U is the 
superposition of all potentials of the nucleons, see eq. (p"5|). 
Let 4> be the potential of a single nucleon in the target. 
The position of the nucleon with number j is denoted by 
the transverse vector rj t ± and the longitudinal coordinate 
Zj. If the range of interaction is much smaller than the 
distance z 1 — z, the potential is practically zero outside the 
domain of integration over £ and the phase shift function 
is given by 



(44) 



X(r x -rj- )X )= / d£(£(r x -T jt ±,^-Zj) 



We have replaced r (£) by r (zj). This means, we use only 
an average value of the transverse coordinate for calculat- 
ing the phase shift for scattering off a single nucleon. 

The two path integrals sum over all possible trajecto- 
ries of the two particles. In order to calculate the cross sec- 
tion for pair production in the nuclear medium, we have 
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to average over all nucleons. We obtain with the path- 
integral (|42| ) 



Gq { r '±,qi z ' l r ^.9i z ) (rj_,gj z' \rx,q, z ) 



Pq-2 . Pq ■ 2 I ,-, 

y r 9 + y x « 1 ^ 



A 

(x (t 5 (z,) - r,- ± ) - X (t, (*,•) - r il± )) | ^ , (45) 

with boundary conditions T q (z) = Tj_,q, T q (z') = r'j_ , 
Tq(z) = r±,q and Tq(z') — r'^ -, The averaging proce- 
dure is similar to the one described in |l9| . We neglect all 
correlations between the nucleons and introduce the av- 
erage nuclear density n A , which is normalized to A, eq. 
([46|). Then, the whole expression may be written as an 
exponential, if A is large enough, 



exp I i *&(z' — Zj)$ (zj — z) 
( X (t 5 (zj) - r jtX ) - X (r q (Zj) - r jtJ 



(46) 



(47) 



= U-\J d 2 s j ' di'n A (s,n (l 

^ Z ^ 

- exp ^ (X (Tq (0 ~ s) - X (Tq ($ ') - s))^ l 

«exp|-y dt'n A (b,0 J d 2 s 

^ Z 

x ^1 - exp (X (Tq (£') - b) - X (Tq (£') - s)) yj I . 

When n A is varying fairly smoothly inside the nucleus, we 
can replace its dependence of s, the transverse distance 
between the pair and the scattering nucleon, by the im- 
pact parameter b, eq. d47j). Since we have a short ranged 
interaction, it is reasonable to choose as impact parameter 
b = (rj_ + r' ± ) /2. This way we find the forward scattering 
amplitude for a dipole scattering off a single nucleon, see 
eq. ([47l). The corresponding cross section, 

*&(p) = 2K J d 2 s fl-exp(i(x(p-s)-X(s))U, 

(48) 

appears as imaginary potential in the propagator. Omit- 
ting the realpart of the amplitude, which is known to be 



small, we finally arrive at 

Gq ( r J_,q> Z ' \ V X : qi z ) G| (?X,q> Z ' l r J-,?j z , 



J Vr q J VT- q exp jz f d^f 



Pq ^.2 . Pq -2 

1 + y T * 



(49) 



+ o«A (biOVqqdTq-Tq 



Note, that all dependence on the potential of the nucleons, 
4> (r), has been absorbed into a qq . Although eq. (|48;) looks 
different for a nonabelian potential, we belief, that the 
formulae presented in this paper still hold, if a qq is taken 
from experiment or calculated in perturbative QCD. 

It is convenient to introduce center of mass coordi- 
nates, T re i = Tq — T q and T cm = (1 — A) T q + Ax g , and to 
express the sum of the two kinetic energies as the sum of 
the kinetic energy of the relative motion and the center of 
mass kinetic energy. 



Gq ( r i,g' z ' l r -L,<ji z ) G\ ( r i,<7' z ' l r -Ugi 2 
2>Tcm / T>T rei exp< i / d£ 



' ff\ i 2 cm 



(50) 



+ 2 UA ( b >^ a qq ( T r<*) 



2ni (z' - z) 



v ( A ( r i, 9 ~ *±«) + (1 - A ) ( r i,<r - r±,q))' 



x exp| i- 



x jVTrei exp{i j d£ ( -T rel + -n A (&,£) a qq (r rel ) 



(51) 



(52) 



We have introduced the reduced mass 
1 _ 1 1 1 

T^P^ + 7q ' ^A(l-A)- 

Since the imaginary potential depends only on the relative 
coordinate, the center of mass propagates freely and what 
remains is the effective propagator 

W(p',z'\p,z) (53) 

= J VT re i exp ji J d£ {^T 2 el - V opt (b, T reh £)) ^ , 

with p' = r' ± y — and p = r±,q — r± tq and the optical 
potential 



V opt (b, p, z) = -^n A (b, z) o q q (p) . 



(54) 
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It fullfills the equation 



dz' 2v\{l-\) 



W(p',z'\p,z) 



iS(z' ~z)5 {2) (p' - p). (55) 



The propagator for the center of mass coordinate produces 
a ^-function and thus, we obtain from eq. ( j39| ) after aver- 
aging over the medium 



7* A 



\v 2 



d 2 b2$t 



dX 



x 2 (i-xy 



dz J dz'e lq ™ n ( z - z ') 

z 

x 0(p)0*(p')W(p , ,z'\p,z) 



p=p=o 



with 



O (p) = O t (p) + O l (p) , 
where the transverse part of the operator is 



(56) 
(57) 



Or (p) = X f q { m f cr ■ e T + i (1 - 2A) er ■ e z e T ■ V (p) 

+ (e, x e T ) ■ V (p)} Xf (58) 

and the longitudinal part 



O l (p) = x\{2 Q X (1 - X)} Xl} = 2 Q A (1 - A) 5 Xq . r 



(59) 

Equation (pq) is the central result of this work. It is 
the total cross section for production of a gg-pair from 
a virtual photon scattering off a nucleus. We have not 
summed over the spins of the quark and the antiquark 
and not averaged over the polarizations of the photon. We 
have not summed over the different flavors, either. The 



expression for the operator 0, eq. (57), depends on the 
spin vector of the quark and the antiquark. The directions 
of these vectors may be fixed arbitrarily. 

In order to represent the result in the form (JXl]) , we 
rewrite W and it's complex conjugate in an expansion. 
The results can be combined in the following way: 

W(p',z'\p,z) 

= W (p',z'\p,z) 

z' 

+ i j ' dzi J d 2 piW a {p',z'\p 1 ,z 1 ) 

Z 

x V op t (b, pi,zi) W (pi,Zx | p, z) 
dzi / d 2 pi I dz 2 I d 2 pi 



x Wo{p\z'\ Pl ,z x ) V: pt (b, Pl , Zl ) (60) 
x W (p 1 ,z 1 \p 2 ,z 2 )V pt{b,p2,z 2 ) W {p 2 ,z 2 \p,z) . 



Here, Wq is the propagator corresponding to (p3|), when 
the potential is absent. The first term gives a divergent 
contribution, which is the wave-function renormalization 
for the photon. The second term leads to the first contribu- 
tion in (p"l|). The operators Ot,l applied to the free propa- 
gator Wo give the light-cone wavefunctions <I>t,l {p, A), up 
to a constant overall factor. The third term in the above 
expansion is the interference term. Since this term con- 
tains the full propagator, there are no higher terms in this 
expansion. 

As an example for further calculation, consider the in- 
tegral 



2 = — 



dz / dz'e tqrn ( z - z ') 



dzi I d 2 pi / dz 2 / d 2 p 2 



x V* t (b,pi,zi) V opt (b,p 2 ,z 2 )W (p',z'\p 1 ,z 1 ) 



x W(p 1: zi\p 2 ,z 2 )W Q (p 2 ,z 2 \p,z), 
that is needed to calculate the interference part 



(61) 



Zja er , 
Av 2 



d 2 b2$ 



dX 



x 2 (i-xy 



o( P )o* (p') i 



p=p=o 



(62) 



With the new variable e 2 = X (1 — A) Q 2 + m 2 f we find 



2u\ (1 - A) 



(63) 



and because of the relation 
8 , A ± (p')-e 2 



dz' 



2vX(l - A) 



(Wo{p',z'\p,z)e-^ n ( z '- z )) 

= i6(z'-z)6W(p'-p), (64) 
we can write the propagator in the form 

W (p\z'\p,z)e-^"' l ( z '- z ) (65) 

OO 

d?l\ f du exp {— icu (z ' — z) + zlj_ • (p f — p)} 



(27T) 



2tt 



2i/A(l - A) 



17] 



The +i?/-prescription for the pole in the complex w-plane 
ensures that Wo (p', z' \ p, z) = for z > z'. Putting ( |6^ ) 
into ( |6l|) yields after a short calculation 

z _ 4^ 2 A 2 (1-A) 2 
(2n) 2 

OO Z\ 

x J dzi J dz 2 e z ^" ,{z2 - Zl) J d 2 Pl J d 2 p 2 

— OO — OO 

x V* pt (b,p 1 ,z 1 ) V opt (b,p 2 ,z 2 ) (66) 
x K Q (e |p' - Pi|) W (p^zi | p 2 , z 2 ) Koie \p 2 - p|) . 
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Kq is the MacDonald function of zeroth order. We have 
used the relation 



photon, A, 



-1 for positive helicity and A 



■1 for 



1 f e il± 'P 



(67) 



We insert (|S6|) into (|62j) and calculate the contribution 
from the second term in the expansion ( |60| ) in a similar 
way. We obtain for the total cross section 



7* A . 7 2 a em 

?t n f — J\ Zj 



/7^T2 J d\ I d pia qq { Pl ) 
o 



O(p)K ( S \p- Pl \) 



p=o 

oo oo 



Jf (2tt) s 



23? / d 2 b / dzi / dz 



-OO Z\ 



dX / d 2 pi / d 2 pi e 



-iq™ m (z 2 -zi) 



X V opt 0, Pi) ^l) Vopt (b, p 2 ,z 2 ) 

x (o*(p')K (e\p'- Pl 
x W(p 2 ,z 2 \p 1 ,z 1 ) 
x (o(p)K (e\p-p 2 \) 

With help of the relation 

V (p) K (ep) = -e^Ki (ep) 
P 

we find the light cone wave functions 



<2> T (sp) = Z f ^^O T (p)K (ep) 



(68) 



(69) 



(70) 



\fa~e 



2tt 



mK Q (ep) Sx^XjSx^ 



+ (iX q (2A - 1) e T • e p + (e T x e z ) ■ e p 
x ^(ep)^, _ A J (71) 



and 



<P L ( £ p) = Z f ^O L (p)K (sp) 



(72) 



= Z 7 ^^2 Q A (1 - A) K (ep) S Xq ,^ (73) 

We made use of the Kronecker-i5. The unit vector in p- 
direction is denoted by e p . As spin vector we have chosen 
the unit vector in z-direction and the parameter X q takes 
the value +1 for spin in positive z-direction and the value 
— 1 otherwise. For the antiquark it is vice versa. For the 



negative helicity. The transverse light cone wave function 
has one part that depends on Kq and another part depen- 
dend on Ki, the MacDonald function of first order. Note 
the different spin structures of these parts. In the Ko-part, 
the spins of the quarks add up to the spin of the photon, 
but in the i^i-part of the transverse light-cone wavefunc- 
tion, the spins of the quarks add to and the pair gets 
a spatial angular momentum. We finally sum over all fla- 
vors, colors, helicities and spin states and get the following 
expression: 



Tf 1 ' A - A 



3a 
(2 



dX f d 2 pa^(p) (\$ T (ep)\ 2 + \$ L (ep)\ 2 



^Y / Z 2 dtjd 2 b J dzx J dz 2 

f ^ — oo zi 



-i<j£"™(z 2 -2l) 



x J dX J d pi J d l p 2 e 



X n A (b, Zi) n A (6, z 2 ) cr^ (p 2 ) cr^ (pi) 
x |(l-2A(l-A)) £ 2 ^-^X 1 ( £ p 1 )^ 1 (ep 2 ) 



PlP2 



4Q 2 A 2 (l 



W(p 2 ,z 2 \p 1 ,zi). 



X) ) K (epi) K (ep 2 ) 



(74) 



Here, \<Pt,l (sp)\ are the absolute squares of the trans- 
verse and the longitudinal light-cone wavefunctions, sum- 
med over all flavors, see eq. (0). This form was used in 
[ p"8| for a calculation of nuclear shadowing. Eq. (Q) was 
for the first time suggested in a paper by Zakharov pq| . 

Let us summarize the assumptions and approximations 
entering this derivation. We start from the Dirac equation 
with an abelian potential. We use this simplification, since 
we know, that all dependence of this potential will be put 
into the dipole cross section. Because the final result inter- 
polates between Aa] ot N and eq. ([)]), we assume, that our 
results hold also for the case of a nonabelian potential. 

Further, no interaction between the quark and the an- 
tiquark is taken into account and therefore, the two Dirac 
equations decouple. We use the Furry-Sommerfeld-Maue 
wavefunctions p2||23[ | that are known to be good approx- 
imations to the continous spectrum of the Dirac equa- 
tion for high energies. Then we derive a two dimensional 
Schrodinger equation for a scalar function that may be 
regarded as an effective wavefunction. The z-coordinate 
plays the role of time, since the particles move almost with 
the velocity of light. This Schrodinger equation is solved 
in terms of it's Green function. Averaging over all scatter- 
ing centers in the nucleus yields an optical potential that 
is proportional to the total cross section for scattering a 
qq-pair off a nucleus. We have omitted the real part of 
the forward scattering amplitude. All dependence on the 
potential is absorbed into this cross section. We propose 
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to use the cross section as input for calculations and not 
the potential from which it originates. It may be taken 
from experimental data and is the nonperturbative input 
for our formulae. Analysis of hadronic cross sections p6| 
suggests <Jqq{p) ~ Cp 2 with C between 2.5 and 3. Note 
however, that we do not make any assumptions on the 
shape of cr^(p) in the derivation. 

The averaging procedure and the summation of the 
multiple scattering series is similar to the one in Glauber 
theory [l9| ] and most of the approximations come in at this 
point. First, we neglect all correlations between the nucle- 
ons. Then, the influence of the potential is described by a 
phase shift function. We also assume, that we have a short 
ranged interaction and the pair interacts only with one nu- 
clcon at a given time. It has been demonstrated by Mueller 
ifjof . that in Born approximation the dominant contribu- 
tion comes from graphs, where the qq-pair interacts with 
the different nucleons one after another via two gluon ex- 
change. Graphs with crossed gluon lines are suppressed. 
This observation justifies our summation procedure. The 
phase shift for scattering a particle in the pair off a single 
nucleon is calculated for an average value of the transverse 
coordinate of the particle. This means, the transverse co- 
ordinates should not vary too rapidly within a longitudinal 
distance of the order of the interaction range. When we 
want to obtain an exponential from the averaging proce- 
dure, the nuclear mass number A should be large enough. 
Further approximations are, that both particles in the pair 
see the same nuclear density. We use the value of the den- 
sity in the middle between the quark and the antiquark 
for our calculation. We also approximate the motion of 
the center of mass of the pair by a free motion, since the 
pair is scattered predominantly in forward direction. 

We finally arrive at the result eq. (^6|). This formula al- 
lows to calculate the cross section a^ ot for arbitrary polar- 
ization of the photon and the pair. However, this equation 
is not convenient for numerical calculations and we mod- 
ify our result, introducing the light-cone wavefunctions, 
eq. (|6§|). Since we are not interested in certain polariza- 
tions for a calculation of nuclear shadowing, we sum over 
all helicity and spin states, arriving at eq. (|74|). 



not important. In particular, the structure function F2 of 
the proton, calculated in our model, does not depend on 
xb ■ In order to get the steep rise at very small x b , one has 
to take higher Fock-states of the photon into account. Fur- 
ther, shadowing for the longitudinal cross section drops as 
1/Q 2 as Q 2 —* 00, since we have no gluon shadowing in 
our model. Adding one gluon to the qq-pair would change 
this. This problem will be addressed in a forthcoming pa- 
per. 

All assumptions summarized in the end of sec. || seem 
reasonable to us and the approximations should work as 
long as the nuclear mass number A is not too small. 

Numerical calculations |l8| with eq. ( f74| ) show that 
higher order scattering terms have a significant influence 
on the total cross section, especially for heavy nuclei. In 
[ p"8[ , the formula for nuclear shadowing was given without 
derivation. This has now been made up. 

The cross section cr^ ot A is identical to the total cross 
section for production of a qq-pair from the virtual pho- 
ton in the field of the nucleus. The suppression occurs 
because of destructive interference between pairs created 
and different longitudinal coordinates within the coher- 
ence length. Thus, shadowing may be regarded as the 
Landau-Pomeranchuk-Migdal-effect |^7],^8| for pair pro- 
duction. This effect is the analog of the more widely known 
effect for brcmsstrahlung and was first mentioned by Mig- 
dal (28) for electron-positron pair production in condensed 
matter. However, this effect is practically not observable, 
because of the low density of solids. Since the density of 
nuclear matter is much higher, the effect occurs in qq pair 
production in DIS. For this reason, eq. (|74|) was discov- 
ered independently by Zakharov p5| , who considered the 
LPM-effect for finite size targets. 

We are grateful for stimulating discussions to Jorg Hufher and 
Boris Kopeliovich who read the paper and made many useful 
comments. We thank the MPI fur Kernphysik, Heidelberg, for 
hospitality. 
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3 Summary and conclusions 

We have considered nuclear shadowing in the rest frame of 
the nucleus, in which the virtual photon fluctuates into a 
qq-pair. In the preceding section we gave a detailed deriva- 
tion of a formula for nuclear shadowing in DIS that ac- 
counts for both, the finite lifetime of the hadronic fluctua- 
tion and all multiple scattering terms. With this formula, 
eq. (|74|), it is possible to calculate nuclear shadowing for 
moderate values of xb, %b > 0.01, where the lifetime of 
the fluctuation does not exceed the nuclear radius by or- 
ders of magnitude. It must be clearly emphasized, that 
we have only taken the qq-Fock component of the photon 
into account and thus, the applicability of our results is 
restricted to values of xb and Q 2 , where corrections from 
higher Fock-states of the photon, containing gluons, are 
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